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SOME  PROBLEMS  IN  THE  DESIGN  OP  THE  TUBE  WIND  TUNNEL* 

Wang  Sung-gao 
(Institute  of  Mechanics) 

ABSTRACT 

The  problem  of  the  boundary  layer  growth  in  the 
charge  tube  is  discussed.  An  analytical  solution  is 
derived  and  has  been  reduced  into  an  algebraical  express 
ion.  The  result  contains  various  factors  which  affect 
boundary  layer  growth.  The  calculation  is  simple  and 
the  results  agree  with  experiments. 

Based  on  P.  L.  Shope’s  model, analytical  expressions 
for  the  test  section  starting  process  are  derived. 
Various  factors  which  affect  the  starting  process  are 
analyzed.  The  precision  of  calculation  is  adequate 
for  design  purposes.  The  calculating  procedure  is  much 
simpler  than  F.  L.  Shope’s. 


NOTATIONS 
a  speed  of  sound 
d  diameter  of  tube 
L  length 
M  Mach  number 
p  pressure 
r  radius 
T  temperature 
t  time 

u  axial  component  of  velocity 

V  velocity  of  expansion  wave 
or  volume 

x  axial  coordinate 

y  radial  coordinate 

Cf  coefficient  of  surface 
friction 

f - 

received  in  October  1981. 

1 


Re  Reynolds '  number 

Y  adiabatic  index  for  gases 

t  characteristic  time  of  open¬ 

ing  of  ^valve 

A  area  or  coefficient 

Subscripts 

3  parameter  of  central  flow  in 
gas  storage  tube 

4  state  of  stored  gas 

e  outward  flow  parameter  of 
boundary  layer 

t  parameter  of  the  test  section 

w  wall  parameter 

p  storage  chamber  parameter 

pe  exhaust  parameter  of  storage 
chamber 


5  thickness  of  boundary  layer  f  parameter  of  regulating  plate 

6*  displacement  thickness  of  s  state  after  starting  time 

boundary  layer  eO  parameter  at  the  stationary 

0  momentum  thickness  of  point  of  the  flow 

boundary  layer 

6f  width  of  opening  of  the 
r  regulating  plate 
0  density 

The  tube  wind  tunnel  has  a  simple  structure  and  it  is  relative- 

* 

ly  easy  to  increase  the  pressure  of  the  stored  gas.  As  the  ratio 
of  the  pressure  at  the  stationary  point  of  the  flow  in  the  test 
section  to  that  of  the  stored  gas  is  determined  by  the  ratio  of 
the  area  of  the  stored  gas  tube  to  that  of  the  nozzle  throat,  one 
can  increase  the  Reynolds'  number  of  the  flow  in  the  test  section 
by  raising  the  pressure  of  the  stored  gas.  Hence,  sufficiently 
high  Reynolds '  numbers  can  be  obtained  with  small  equipment .  The 
Reynolds'  number  can  be  further  increased  by  combining  this  method 
with  low  temperature  techniques. 

Some  special  problems  are  associated  with  design  of  tube  wind 
tunnels  to  be  used  as  high-Reynolds '  number  transonic  test  equip¬ 
ment.  We  will  discuss  two  of  these  problems  in  what  follows: 


1.  Boundary  layer  growth  in  the  gas  storage  tube 


In  the  tube  wind  tunnel,  a  one-dimensional  stable  constant  flow 
is  created  during  the  process  of  propagation  of  unstable  non¬ 
constant  expansion  waves.  Owing  to  the  viscosity  of  the  flow,  a 
boundary  layer  is  created  along  the  tube  wall,  its  thickness 
increasing  with  time.  This  boundary  layer  growth  affects  the 
effective  flow  duration  and  quality  of  flow  in  the  tunnel. 

This  problem  has  been  discussed  by  E.  Becker  [1],  J.  C.  Sivells 
[2]  and  H.  Ludwieg  [3]-  All  of  them  used  complicated  numerical 
methods.  J.  C.  Sivells  had  to  modify  his  results  on  the  basis  of 
experimental  results.  Our  analysis  consists  of  directly  integrating 
the  integral  equation  of  the  momentum  of  the  unstable  non-constant 


boundary  layer  and  reducing  the  analytical  solution  thus  obtained 
to  an  algebraic  expression.  The  axially  symmetric  integral  equa¬ 
tion  of  momentum  takes  the  following  form  for  cylindrical  flows: 

1  9  .  \  99  _  1  »  Cl) 


where 


1.  ~  —  Ct 

dt  dx  2 


-I  ('-*)£('—)* 
*•- KC'-^'-^)* 


Let  §  —  Vt  —  * 

Here,  V  is  the  velocity  of  the  unstable  non-constant  expansion 
waves.  Its  relation  with  the  speed  of  sound  a^  in  the  stored  gas 
state  has  already  been  derived  by  E.  Becker: 

„  4-(r  +  l)M» 

—  rr  +  rr~  cm 

«•  1  +  M, 

2 

is  the  Mach  number  corresponding  to  ug. 


Assume  that 


?*J=A 

e 


is  independent  of  time,  and  we  have 


V  «.  9  /  de  2  ' 


\  •,  9  l  JS  2 

We  use  equation  (5)  to  analyze  the  following  factors: 

(1)  Assuming  that  Cf  can  be  expressed  as  a  power  expansion  in 
0,  and  that  the  velocity  distribution  can  be  expressed  as  a  power 
expansion  in  we  can  derive  from  » >  7,  and  */?$  <i  the  relation 

«•/•*  <  10/7 

Here,  n  is  the  inverse  of  the  index  of  the  expansion  of  the  velocity 
section. 


(2)  The  effect  of  Cj.  is  obvious.  Different  expressions  have 
been  given  for  C-  for  different  Reynolds'  numbers  [4].  Comparing 


the  formulas  given  by  Blasius  and  Karman,  the  result  obtained  from 
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one  is  more  than  twice  that  obtained  from  the  other,  for  Re  >  10'. 

(3)  The  effect  of  the  width  of  the  expansion  waves  is  on  the 
velocity  V  of  the  effective  expansion  waves.  It  can  be  seen  from 
equation  (4)  that  this  effect  can  be  neglected  only  for  small 

(4)  The  effect  of  the  velocity  section  is  expressed  through 
n.  Taking  a  two-dimensional  slab  as  an  example,  we  have  approxi¬ 
mately 


6_ _ Z - , 

8  («  +  1)0»  +  2) 

n  varies  as  Re  [51. 


In  the  process  of  solving  equation  (5),  we  have  been  able  to 
take  all  of  the  above  four  factors  into  account.  Calculations  show 
that  when  n  *  7-9$  the  error  pro.duced  by  the  assumption  that 
(**  — *,)/*  is  independent  of  time  is  approximately  which  is 
acceptable . 


We  take  the  following  expression  for  the  coefficient  of  fric¬ 
tion.  This  is  a  combination  of  Karman' s  and  Frankie,  F.  Voishel, 
V’s  (sic)  equations. 


C, 


- - -(l  +  Z^iMlV 

(lgRe«  +  1.1696)(lg Re#  +  0.3010)  '  2  / 


(6) 


Inserting  this  in  equation  (5),  and  integrating,  we  obtain, 
after  rearrangement, 


After  — ,  Can*#* 


ff_  0.0293 {Vt-x') 

e  +  v 

9  u.  8 


flg(2Re#)J“,^I  +—-J-MJ) 


(7) 


are  calculated,  the  above  equation  can  be  written 
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as 


C  -  B(«/r.)  -  A(8/rty 


(8 


where 


A  — 


_ (i+-L-Il)  +  I - ” - 

+  3)\  2a  +  1  TJ  u.  (2n  +  l)(2a  +  2)  V  »  Tj 


(2a  +  2)(2a 

b _ 2 _ (h-L-LU  L _ " _ 

(a  +  2)(a  +  3)  \  a  +  1  Tm>  u,  (a  +  l)(a  +  2)  V  »  Tm / 

c  _ _ 0.0293(r/  -  ») _ 

r.tlg(2R*)]’(l  +  -2L=LLMi)*«  c- 


(9 


When  we  are  calculating  for  the  thickness  of  the  boundary 
layer  at  the  entrance  of  the  nozzle,  x  =  0.  If  we  let  ■ 
then 

- • - +  L _ i _ 

(2a  +  l)(2a  +  2)  u,  2(2 a  +  1) 

»  y  V  1 

(a  +  l)(a  +  2)  u,  a  +  1 


B 


h|h 


The  procedure  for  the  calculation  is  very  simple:  given  Tq, 

M-,  Re.  and  a,..  Assume  values  for  —  and  find  Re*  and  Re..  Prom 

3  d  v  4  vr,  60 

find  —  and  V.  From  n  and  — ,  find  A  and  B.  After  C  is  found, 
canebe  readily  obtained.  eFigures  1  and  2  give  the  calculated 

rt 

results  as  compared  with  experimental  results.  The  agreement  is 


very  good.  Our  computation  has  a 

for  i-<0.75 
u 


Figure  1.  Comparison  of  cal¬ 
culated  results  with  experimental 
results  of  AEDC 
+  experimental  values  for 
M3*0.265  Re  *3 • 0x10 7 

- J.  C.  Sivells  calculation 

- our  calculation 


fairly  high  accuracy,  at  least. 


Figure  2.  Comparison  of  calculated 
results  with  experimental 
results  of  Goettingen 
H  experimental  values  for 
M3*0.2  Red  5 . lxlO4 

o  experimental  values  for 
M3*0 . 3  Red=7.3xl04 

- H.  Ludwieg’s  calculation 

- our  calculation 


2 .  Analysis  of  the  starting  process 

The  tube  tunnel  is  an  equipment  for  short-duration  operations. 
The  total  process  takes  place  within  a  time  period  on  the  order  of 
10  msec,  while  the  starting  time  for  the  flow  in  the  transonic  test 
section  to  reach  equilibrium  can  sometimes  be  on  the  same  order  of 
magnitude,  and  takes  up  a  sizeable  portion  of  the  total  time.  To 
increase  the  effective  operation  time,  one  must  shorten  the  starting 
time.  Reducing  the  starting  time,  however,  results  in  an  increase 


in  the  aerodynamic  load  on  the  wall  of  the  test  section  and  the 
regulating  plate.  Therefore,  it  is  important  to  study  the  starting 
process,  analyze  the  major  factors  and  obtain  simple  and  convenient 
computation  methods. 


We  employ  the  method  of  concentrating  parameters  [6].  We 
found  out  from  analysis  that  an  analytical  solution  can  be  obtained 
by  using  an  ordinary  differential  equation  to  describe  *  starting 
process.  Under  the  conditions  of  supersonic  flow,  the  nation  of 
conservation  of  mass  within  the  storage  chamber  is 


-i-  ( ItS]  4.  Ah  (  ?*  \  At  —  A, 

dt\pH)  1.73  V,  \pM/  1.73  V , 


The  solution  satisfying  the  initial  conditions  is 

pa  ' 


(n: 


(12: 


At  the  end  of  the.  starting  process,  p 


P 


Alpt* 


and  we  obtain 


_£a_  Ah.  f*  —  in  +  (At — 4»VAtr, 

1-73  V,  Ajt,  +  ( A #  —  A,)/ Ah 


(13: 


In  the  above  equation, *-<*/*)  U,  A,  reflect  the 
resistive  effect  of  the  wall. 

M  ^1  &£%  — _  4 \ 

When  Ap  *  0,  we  have  pM  ***  ah  1.73  v  ’  •  Except  for 

A  *  A^,  there  is  no  balance  between  the  two  sides  of  the  above 
©  t 

equation. 


For  transonic  flows,  there  is  no  closed  ordinary  differential 
equation  that  can  be  used.  On  the  basis  of  experimental  and  cal¬ 
culated  results,  we  assume  that  [73 


Al.  —  ^  4.  (1  — 
f* 

*-U4t-(t  +  f)  | 


Let 


The  mass  conservation  equation  in  the  storage  chamber  is 


Its  solution  is 


±.  (-&-)  +  N  II.  -  K'  [?  +  (1  - 

at  \  pA  /  1.73  V , 


It. 

Pi 


K' 


1  +  K 


K ' 

1  +  K 


_ d_  > 

A ft. 


(16) 
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K' 


1  +K 


1 


K' 


1  - 


e  i 

Nr. 


(17) 


When  t  =  ts,p  *  ^ipt»  and  e9uation  (17)  is  a  transcendental  equa¬ 
tion  which  can  only  be  solved  by  numerical  methods.  However,  for 
the  case  where  r,-*o  and  the  thin  plate  is  used  for  starting,  or 
for  the  case  where  r.-*oo  and  the  valve  is  used  for  starting,  we  can 
obtain,  respectively 


and 


(1+K) 


af  Ap, 
1.73  V, 


I,  1  — JK'/Q  +  K)]ft 
ft[Al~K70  +K)] 


(18) 


0  +10 


ar*  Ap. 

1.73  V, 


U 


i  -  [xyq  +k)] 
Au*  -  K'/O  +  K) 


(19) 


To  obtain  meaningful  results,  we  must  place  restrictions  o.n 

KV(1+K).  This  places  restrictions  on  the  relations  among  A^Ay, 

and  A  .  In  other  words,  the  geometrical  parameters  of  the  storage 
Pe 

chamber  are  interrelated. 


Similar  equations  and  relations  hr. Id  for  the  case  where  Ape  =  0. 
Analysis  shows  that  only  when  A^.  is  large  and  there  is  sufficiently 
strong  induced  emission  in  the  pressure  expansion  section  will  the 
equipment  start  smoothly. 

Prom  the  above  analysis,  we  know  that  the  major  factors  affect¬ 
ing  the  starting  process  ar eA*/A„  Vf/V„  A,/A,*  A./A,  and  Lt/a..  .  With 
respect  to  the  tube  wind  tunnel,  the  auxiliary  exhaust  is  very 
important,  and  the  duration  and  frequency  of  opening  the  exhaust 
has  a  very  large  effect  on  the  starting  load.  A,  —  a,  c  Ap.  is  a  con¬ 
dition  favorable  for  the  flow. 
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We  performed  calculations  on  several  published  sets  of  data 
for  AEDC  model  tube  wind  tunnels  using  our  analytical  method  of 
solution.  The  method  is  easy  and  convenient  to  apply  and  except  for 
certain  special  cases,  the  results  obtained  agree  very  well  with 
those  obtained  via  numerical  methods.  The  discrepancies  are  gener¬ 
ally  10%  or  less.  Por  example,  for  the  case  of  starting  with  the 
value  where  M  *  1.10,  the  AEDC  value  t_  =  76  ms,  while  the  analy- 
tical  solution  gives  t_  =  80  ms. 

S 
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SIMPLIFIED  NAVIER-STOKES  EQUATIONS  AND  COMBINED  SOLUTION 
OF  NON-VISCOUS  AND  BOUNDARY  LAYER  EQUATIONS* 


Kao  Chi 

(Institute  of  Mechanics,  Academia  Sinica) 

ABSTRACT 

This  paper  presents  a  part  of  the  technical  report 
[2]  in  which  the  author  studied  the  simplified  Navier- 
Stokes  equations  and  the  combined  solution  of  nonviscous 
and  boundary  layer  equations.  From  the  full  Navier- 
Stokes  equations  and  an  analysis  of  the  combined  solution 
of  nonviscous  and  boundary  layer  equations,  simplified 
Navier-Stokes  equations  were  worked  out .  A  perturbation 
analysis  which  differs  slightly  from  the  match-perturba¬ 
tion  expansions  of  inner-outer  layers  developed  by  Van 
Dyke  [1]  shows  that  the  solution  of  the  simplified  Navier- 
Stokes  equations  is  uniformly  valid  with  accuracy  of 
0(Re  in  the  whole  flow  field,  where  Re  =P«>-U«>L,  p  * 

w  oo  y  "  oo 

is  the  density  of  free  stream,  U^  the  x-componen£  of  velo¬ 
city,  L  the  characteristic  length,  yo  the  dynamic  viscosity 
of  free  stream. 

The  simplified  N-S  equations  possess  the  properties  of  para¬ 
bolic-hyperbolic  equations.  Under  equilibrium  conditions,  it  is  much 
easier  to  use  the  method  of  the  forward-progressing  calculation  to 
solve  the  simplified  N-S  equations  than  to  solve  the  full  elliptical 
N-S  equations  by  means  of  numerical  methods.  .While  solving  the 
simplified  N-S  equations,  one  simultaneously  obtains  the  non-viscous 
external  flow  as  well  as  the  viscous  boundary  layer  flow.  Theoret¬ 
ically,  this  is  superior  to  the  conventional  procedure  of  first  com¬ 
puting  for  the  non-viscous  flow  and  then  computing  for  the  viscous 
boundary  layer  flow.  With  respect  to  many  types  of  flow  fields,  the 
simplified  N-S  equations  can  realistically  reflect  their  mechanical 
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aspects.  For  instance,  one  can  use  these  equations  to  accurately- 
calculate  for  the  complex  flow  diagram  of  the  interaction  among 
viscous  boundary  layer,  the  high-entropy  layer  and  the  non-viscous 
external  flow  of  hypersonic  diffracted  flow  fields. 

Methods  for  solving  the  simplified  N-S  equations  are  currently 
under  research.  In  fact,  satisfactory  answers  have  not  yet  been 
found  for  problems  regarding  the  mathematical  properties  and  stabi¬ 
lity  of  these  equations,  the  correct  way  of  presenting  the  Cauchy 
problem,  and  the  degree  of  accuracy  that  can  be  reached  b  y  the  flow 
calculations.  Besides,  different  authors  hold  different  points  of 
view.  This  paper  presents  part  of  the  technical  report  [2]  the  author 
gave  in  1967  in  which  he  obtained  the  simplified  N-S  equations  from 
an  analysis  of  the  combined  solution  of  non-viscous  flow  and  viscous 
boundary  layer  equations,  and  used  the  perturbation  method  to  show 

that  one  can  obtain  from  that  set  of  equations  a  solution  that  is 

-1/2 

accurate  to  the  order  CKRe^  ).  It  should  be  pointed  out  here  that 
the  author  later  became  aware  of  a  similar  derivation  of  the  simpli¬ 
fied  equations,  given  in  [3]. 


1.  Simplified  N-S  equations 

Let  x  and  y  be  the  orthogonal  coordinates  along  the  wall  and 
perpendicular  to  the  wall,  respectively  (see  Figure  1).  Let  u  and 
v  be  the  corresponding  velocity  components  and  p,  p,  T,  y  and  A  be 
the  density,  the  pressure,  the  temperature,  the  viscosity  coeffi¬ 
cient  and  the  coefficient  of  thermal  conduction,  respectively.  Take 

2 

the  parameters  and  pw  of  the  oncoming  stream  to  be  the  char¬ 
acteristic  values  of  velocity  and  pressure,  respectively,  and  the 
fixed  wall  length  L  to  be  the  characteristic  value  of  length.  Assume 
that  the  order  of  magnitude  of  p  and  y  can  be  estimated  from  pa 
and  y^.  From  the  standpoint  of  non-viscous  external  flow  and  that 
of  viscous  boundary  layer  flow,  we  make  an  estimate  of  the  various 


terms  in  the  full  N-S  equations  as  follows 
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Here  *  — +  A,  y  —  o  (two-dimensional  flow),  1  (axially  symm¬ 
etric  flow);  H  =  1  +  ky,  where  k  is  the  curvature  of  the  wall  sur¬ 
face;  Re-  —  -Bm,VmL  t  A  _  0(Rezl/r)  where  6  is  the  thickness  of  the 

Mm  L 

boundary  layer;  $>.  and  $  represent  all  the  terms  in  the  tangential 

u  n 

and  normal  momentum  equations,  respectively,  that  are  of  order  of 
magnitude  equal  to  or  less  than  n  J5y  .  We  discuss  the  case  for 
Re->1  as  follows:  L 


1)  For  y  >6,  omitting  the  terms  in  the  full  N-S  equations 

that  are  of  order  o^Rei1-^  we  obtain  the  Euler  equations.  For 
y<8,  ,  we  omit  the  terms  in  the  full  N-S  equations  that  are  of 

order  equal  to  or  less  than  o^R*-'*^— )  and  obtain  Prandtl's  boundary 
layer  equations. 

2)  When  solving  the  Euler’s  equations  and  the  boundary  layer 

equations  simultaneously,  a  mathematical  singularity  exists  on  the 
boundary  line  y  ■  6.  The  asymptotic  value  of  the  solution  of  the 
Euler’s  equations  as  y(>5)  <5  will  not  agree  completely  with  that 

of  the  solution  of  the  boundary  layer  equations  as  y(<6)  6.  The 

actual  order  of  magnitude  of  this  singularity  is  o(KtZM-^j  .  One 
should  note  that  as  the  various  inertial  terms  in  the  normal 

Euler  equation  are  also  of  the  order  of  magnitude 

3)  Therefore,  in  order  to  eliminate  the  singularity  that  arises 

in  the  simultaneous  solution  of  the  non-viscous  external  flow  and 
the  viscous  boundary  layer  flow,  we  should  omit  all  terms  in  they<j, 
full  N-S  equations  that  are  of  order  equal  to  or  less  than  o^Re7‘-^) 
The  resulting  equations  are  the  simplified  N-S  equations.  Solving 
these  equations  simultaneously  with  the  Euler  equations,  we  can 
remove  the  singularity  and  obtain  solutions  with  accuracies  on  the 
order  of  that  are  uniformly  valid  throughout  the  entire 
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flow  field. 


As  described  above,  the  simplified  N-S  equations  are  obtained 
by  neglecting  all  terms  in  the  full  N-S  equations  that  are  of  order 
equal  to  or  less  than  O^Re:*-^2)  as  estimated  from  the  standpoint 
of  the  boundary  layer.  These  equations  are  given  below  (along 
with  the  continuity  and  energy  equations). 


where 


(1.3) 


(1.4) 
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(1.5) 


(1.6) 


c  is  the  specific  heat  under  constant  pressure.  The  simplified  N-S 

r  a] 

equations  (1.3)“ (1.6)  do  not  contain  ~  terms.  When  y  +  0,  they 

Or 

have  four  redundant  characteristic  values  and  two  non-zero  real 
characteristic  values.  Therefore,  these  equations  possess  the  char¬ 
acteristics  of  parabolic-hyperbolic  equations  [4]  and  the  method  of 
forward-advancing  calculation  can  be  applied  under  equilibrium  con¬ 
ditions.  It  should  also  be  noted  that  in  the  simplified  N-S  equa¬ 
tions,  when  y(<6)  6,  ^  from  O(IUJP)-* o(l)  and  the  order  of 

magnitude  of  all  viscous  terms  become  q  .Saj.  Ther^?^’e»  the 
simplified  N-S  equations,  accurate  to  the  order  °(Re*1/*“jif)  »  will 
smoothly  change  into  a  set  of  Euler  equations  for  .  In  other 

words,  to  obtain  the  solutions  for  the  entire  field  that  are  accurate 

13 
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to  the  order  ,  it  is  not  necessary  to  SQlve  the  simpli- 

fled  N-S  equations  simultaneously  with  the  Euler's  equations.  In 
principle,  it  will  be  sufficient  to  solve  only  the  simplified  N-S 
equations,  given  the  initial  and  boundary  values,  i.e.,  the  proper¬ 
ties  of  the  oncoming  stream  (or  shockwave)  and  the  wall  surface. 


2 .  Discussion  on  accuracy 


A  uniformly  valid  and  consistent  solution  on  the  order  of 
0  can  be  obtained  when  one  applies  the  simplified  N-S 


equations  to  the  computation  of  the  entire  flow  field,  including  the 
non-viscous  external  flow  and  the  viscous  boundary  layer.  In  order 
to  prove  the  validity  of  this  statement,  we  need  to  make  use  of  a 
perturbation  expansion  that  is  slightly  different  from  the  usual 
match-perturbation-expansions  of  the  inner-outer  layers  [1].  Within 
the  viscous  boundary  layer  o<y<4  ,  the  solution  of  the  full  N-S 

equations  (1.1)  and  (1.2)  can  be  expanded  in  terms  of  a  — R*:*'1  as 
follows : 

m/Um  *■“#«■+•  itta  +  8*ii„ 

*/ tUm  ™  +  twu  +  «*»#  +  •  •  • 

p/p~  —  p*  +  *Pn  +  *2Pn  ■+■••• 
p!  Pm  —  Pi,  +  a  Pi,  +  a2pit  +  •  •  • 

”  Pm(X)]/apmU a  *■?(»  +  tpn  +  1%  +  ••• 


(2.1) 


in  the  above  X  —  x/L ,  Y  —  y/aL ,  a  —  Rt?'*,  —  —  0(1); 

OY 


and  the  quantities  uio,uii**'*  and  their  partial  derivatives  with 
respect  to  X  and  Y  are  all  of  order  0(1).  Substituting  the  expan¬ 
sions  in  equation  (2.1)  into  the  N-S  equations  (1.1)  and  (1.2)  and 
comparing  terms  in  the  same  power  of  e,  we  obtain 


<><••>  «.(..*?• 


+  r„  -  -  'flrW  +  JL  (^2*0.) 

Ma  - 

+  +  +  wJ  +  pjp* 


(2.2) 


'  dX  dX  )  dX  dY  V  av  +  *“ 


dY 


dY 


) 


P>*  (*'•  qx  +  r"  ^dY  +  ktY*£}  — 

+  ±  JL +  _»_ /  a«.« \  ,  -tiii a  du,t 

3  8Y  V*  dY  )  +  dx  \**  dY  )  +  ;~7~  **Hy 


(2.3) 


In  the  above,  — $l,  r,  —  r/L,  ?.(X)  —  pm(X)/pmUl;  .  Similar  expansions 
can  be  obtained  for  the  continuity  equation  and  the  energy  equation. 
It  should  be  pointed  out  that  in  the  present  expansion,  the  first- 
order  equations  (2.2)  of  order  0(1)  and  the  second-order  equations 
(2.3)  of  order  0(e)  of  the  N-S  equations  for  are  different  from 

those  obtained  from  the  usual  match-perturbation-expansions  of  inner 
outer  layers  [1].  In  the  usual  expansions,  the  first  and  second 
order  equations  for  the  normal  momentum  for*<3  are,  respectively: 


—  0 
dY 


*tiL- 

dY 


Hence,  in  the  usual  expansions,  a  mathematical  singularity  of  order 
O  exits  on  the  boundary  y  *  6  when  either  the  first  order  equations 

or  the  second  order  equations  of  the  inner  and  outer  layers  are 
solved  simultaneously. 


In  the  non-viscous  external  flow  region  where  y  >  5,  the  solu¬ 
tion  of  the  N-S  equations  (1.1)  and  (1.2)  can  be  expanded  in  terms 
of  e  as  follows: 


m/Um 

*/Um 

p/Pm 

fix,  Y,)/pmUi 


■<•  +  ««*  +  t*um  +  •  •  • 

+  •  •  • 

Pm  +  *Pm  +  t*Pm  +•••• 
Pm  +  *pm  +  **Pm  +  ••• 


(2.4) 


In  the  above  equations  Ye-y/L,  and  the  quantities  ueQ,  uel>...  and 
their  corresponding  partial  derivatives  with  respect  to  X  and  Ye 
are  all  of  order  0(1).  Substituting  (2.4)  into  equations  (1.1)  and 
(1.2),  and  comparing  terms  in  the  same  power  of  e,  we  obtain 
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H  ax 


H  QX 


+  —  («#!*'*  + 


H  dx  H  9X 


+ "  (ff Hi- +  &  ~  -  ~  ^ 

\H  dX  OY.  H  )  6Y. 


(2.6) 


It  can  be  readily  seen  that  in  the  region  y  >  5  of  non-viscous 
flow,  the  0(1)  approximation  of  the  N-S  equations,  i.e.,  the  set 
of  first-order  equations  (2.5),  is  the  Euler  equations,  while  the 
approximation  accurate  to  0(e)  i.e.,  the  second-order  equations 
(2.6),  is  a  modified  solution  that  takes  into  account  the  effect  of 
the  displacement  thickness  of  the  boundary  layer.  Solving  equations 
(2.2!)  and  (2.5)  simultaneously;  one  can  obtain  a  solution  accurate 
to  0(1)  that  is  uniformly  valid  in  the  entire  flow  field.  Solving 
equations  (2.2),  (2.3)  and  (2.5),  (2.6)  simultaneously  via  match- 
perturbation-expansions,  one  can  obtain  a  solution  accurate  to 
0(RmZ,/*)  that  is  uniformly  valid  in  the  entire  flow  field. 


Comparing  the  simplified  N-S  equations  (1.4)  and  (1.5)  as  well 
as  equations  (2.2)  and  (2.3)  with  equations  (2.5)  and  (2.6),  one  can 
prove  that  the  solution  to  the  simplified  N-S  equations  satisfies 
the  following  relation,  taking  u  as  an  example 


y  <«  |~  -  («*  +  •«*)!  <  0(«*  -  W) 


jr  >  #  77-  -  Ou  +  «0  <  (HR.;*) 


(2.7) 


Equation  (2.7)  holds  for  the  other  variables  as  well.  Hence,  the 
solution  to  the  simplified  N-S  equations  is  accurate  to  oCRe;*7*) 
and  is  uniformly  valid  in  the  entire  flow  field. 


3 .  Conclusion 


The  actual  form  taken  by  the  simplified  N-S  equations  will 
vary  slightly,  depending  on  whether  the  viscous  terms  of  order 
0($tZu *)  are  neglected  or  not.  In  general,  such  variations  are 
considered  small  [31-C7].  Nevertheless,  to  obtain  a  solution  that 
is  accurate  to  the  order  0<Re;1/2)  and  uniformly  valid  in  the  entire 
flow  field,  it  would  probably  be  advisable  to  use  the  simplified 
N-S  equations  (1.3)-(1»6)  presented  by  the  author  [2], 

The  simplified  N-S  equations  can  be  used  to  obtain  an  accurate 
description  of  the  complex  interference  phenomenon  among  the  non- 
viscous  flow,  the  high-entropy  layer  and  the  boundary  layer  flow 
in  the  shockwave  layer  [31-C7].  They  can  also  be  applied  to  other 
flow  field  calculations,  such  as  non-equilibrium  flow,  flow  in  a 
tube,  distant  wake  flow,  shock-boundary  layer  interference  and  the 
viscous  diffracted  flow  near  pointed  head  and  blunt  head  three- 
dimensional  bodies-  The  simplified  N-S  equations  have  also  been 
applied  to  the  computation  of  such  complex  flow  fields  as  segregated 
flow,  near  wake  flow  and  compressed  flow  near  corners.  One  should 
note,  however,  that  the  simplified  N-S  equations  possess  the  pro¬ 
perties  of  parabolic-hyperbolic  equations  only  when  the  viscous 
effect  is  dominant,  i.e.,  when  y  /  0.  For  y  >  6  where  one  can 
neglect  the  viscous  effect,  i.e.,  y  *  0,  these  equations  become  the 
Euler's  equations  with  corresponding  characteristic  roots 

■  a1-  C* 

Here  C  Is  the  speed  of  sound  and  M  is  Mach  number.  For  M>1, 
are  real  characteristic  roots  belonging  to  the  hyperbolic  family, 
and  the  simplified  N-S  equations  retain  the  same  mathematical 
classification.  For  M<1,  Xm  are  complex  characteristic  roots, 
belonging  to  the  elliptic  family,  and  there  is  a  change  in  the  math¬ 
ematical  classification  of  the  simplified  N-S  equations.  In  the 
region  where  the  simplified  N-S  equations  are  elliptic  equations, 
the  Cauchy  problem  is  not  applicable.  It  is,  therefore,  necessary 
to  obtain  the  solution  by  an  iterative  method  or  some  other  special 
treatment . 
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